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We use Lee-Suzuki mappings and related techniques to construct effective two-body p-shell inter- 
actions and neutrinoless double-beta operators that exactly reproduce the results of large no-core- 
shell-model calculations of double-beta decay in nuclei with mass number A = 6. We then apply 
the effective operators to the decay of nuclei with A = 7, 8, and 10, again comparing with no-core 
calculations in much larger spaces. The results with the effective two-body operators are generally 
good. In some cases, however, they differ non- negligibly from the full no-core results, suggesting 
that three-body corrections to the decay operator in heavier nuclei may be important. An appli- 
cation of our procedure and related ideas to /p-shell nuclei such as ^^Ge should be feasible within 
coupled-cluster theory. 



I. INTRODUCTION 

Particle physicists hope to learn about the overall 
neutrino-mass scale by observing neutrinoless double- 
beta {OiyfUP) decay [l|. To extract a mass from a lifetime, 
however, one must know the value of the nuclear matrix 
element that governs the decay. For that reason, theo- 
rists have worked hard over the last 20 years to better 
calculate the matrix elements. 

One of the best frameworks for the job at present is 
the nuclear shell model. Good calculations (e.g. Refs. 
0, Q) ^se model spaces of dimension 10'' or larger by 
including the full valence shell. Even these calculations, 
however, omit most of the relevant many-particle Hilbcrt 
space by requiring that most particles remain frozen in 
an inert core and prohibiting even active particles from 
sampling levels above the valence shell. These approxi- 
mations induce error that can in principle be accounted 
for through the use of an effective Hamiltonian and decay 
operator. The literature contains a number of schemes 
for constructing effective operators 0, In practice, 
however, such techniques are generally restricted to the 
effective Hamiltonian, and even there the application is 
through a perturbative scheme whose weaknesses often 
must be supplemented by fitting to spectra. The decay 
operator is usually not corrected at all, except at short 
distances (and sometimes through an overall phenomeno- 
logical multiplication factor). How much are calculated 
matrix-element corrupted as a result? How might one do 
better? 

To begin to answer these questions, we look at the 
QvPP matrix clement in nuclei with mass number A be- 
tween 6 and 10. Such nuclei, of course, do not undergo 
/3/3 decay, but one can calculate the matrix elements 
nonetheless. Moreover, in these nuclei we can carry out 
fairly complete no-core shell-model (NCSM) 0] calcula- 
tions and map their results onto valence-shell (p-shell) 
calculations to construct effective operators that repro- 



duce the full matrix elements exactly. We can then by- 
pass perturbation theory, which is often unreliable (and 
was applied inconclusively to /3/3 decay in Ref. Q), and 
test nonperturbative approximations to the full effective 
operator. Reference [8|, |9| carried out this program for 
charge-conserving electromagnetic transition operators, 
the leading pieces of which are one-body. Here, the 
lowest-order effective decay operator acts on two bodies, 
and we examine the restriction to this leading term. We 
define effective operators that reproduce the exact matrix 
element for the artificial decay of states in ^He to those in 
^Be. Are these operators significantly different from their 
bare counterparts? Can they also reproduce Ov^P matrix 
elements in heavier nuclei, or are three- and higher-body 
effective operators necessary as well? The answers will 
provide a good idea of how much work awaits us in the 
heavier nuclei that actually undergo /3/3 decay. 

In Scc.HIlbelow. we describe the concepts and methods 
we employ. Section IHII presents our results, and Sec. IIVI 
discusses their implications for matrix elements in the 
heavier nuclei that are used in /3/3 experiments. 



II. METHODS 

In the closure approximation (which is good for neutri- 
noless decay) and the usual assumption that the nuclear 
weak current is adequately represented by a one-body op- 
erator, the 01^/3/3 matrix element is a sum of three terms: 

Mf. ^ (/I Mg^ + AC + Mj, I*) , (1) 

ab 

the last of which is a very small tensor piece [l^ that will 
be ignored here. The other two Af 's arc given by [Tl| 

M^^ =HGT{rab)<Ta-<Tb (2) 
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with the labels a and b indicating nucleons both here and 
in Eq. ([1]), Tab representing internucleon distance, and the 
"neutrino potentials" H defined by 

HK(r) = — / — dq, K^GT,F. (3) 

nr Jo q + u} 

The hpc (q) in Eq. ([3]) contain the vector and axial- vector 
coupling constants, form factors that account for the 
finite size of the nucleon, and the effects of forbidden 
currents (weak magnetism and the induced pseudoscalar 
term) . The quantity Q is an average intermediate-nucleus 
excitation energy to which the Hk are not very sensitive. 
The authors of Ref. [l^ recently applied chiral effective- 
field theory to derive two-body corrections to the weak 
current and thus three-body corrections to the operators 
in Eq. ([2]). From our point of view, these corrections 
modify the bare 0iy/3(3 operator and are subject to the 
same nuclear-structure renormalization that we apply to 
the two-body operators in Eq. We neglect the chiral 
corrections here to keep matters simple. 

The matrix element Ai fi is often small because of can- 
cellations among contributions at different internucleon 
distances Tab- Instead of looking only at the matrix ele- 
ment, therefore, we also examine the internucleon matrix- 
element distribution C(r), defined e.g. in Ref. [l^l, so 
that 



C(r) dr^Mf^. 



(4) 



Our starting point for calculating matrix elements is 
the NCSM. We use different starting interactions — the 
CD Bonn potential Q and the N^LO chiral effective- 
field-theory interaction [H] — and model spaces that 
allow between six and ten Hlo of excitation energy out- 
side the p shell (i.e. the NCSM parameter N^ax is be- 
tween 6 and 10). We first apply standard Lee-Suzuki 
techniques p^l - ll8l | to the Bonn potential and the Simi- 
larity Renormalization Group (SRG) [l3,[23] to the chiral 
potential to construct interactions appropriate for those 
model spaces. In principle the double-beta decay op- 
erator should be treated in the same way. Preliminary 
studies 
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show, however, that the renormalization is 
slight and confined to short distances as expected, and 
instead of carrying it out here we simulate short-range 
effects through an effective Jastrow function from Ref. 

Many of the isotopes we discuss are very weakly bound 
or unbound in reality, and our representation of them in 
the oscillator-based NCSM distorts their structure. For 
our purposes, however, the poor representation is not 
important; we want to examine the effect of moving to a 
much smaller model space, and take the large-space cal- 
culations to be the "exact" results we want to reproduce. 

Our small model space consists of all but four particles 
residing anywhere in the Op shell and the rest forming an 
inert Os-shell core. As in Refs. [1, 0], we first equate 
the effective neutron P3/2 and single-particle ener- 
gies to the two lowest-energy eigenvalues produced by 



the full calculation in ^He, and the effective proton ener- 
gies to the corresponding eigenvalues in ^Li. Then in the 
A ~ 6 nuclei we use the Lee-Suzuki procedure to map 
the two lowest J'^ = 0+ states, the lowest l"*" state, and 
the two lowest 2+ states (all with T = 1) onto corre- 
sponding orthogonal p-shell states. In doing so we have 
assumed isospin conservation in our small-space calcula- 
tion; breaking isospin would require only the additional 
straightforward step of carrying out separate calculations 
in He and Be. 

The Lee-Suzuki mapping, which comes as close as pos- 
sible to making the p-shell energy eigenstates the projec- 
tions of the corresponding full-space states without spoil- 
ing orthogonality, proceeds as follows. We let P project 
onto the c?-dimensional small space, let Q = 1 — P, and 
denote by \p), \p'), \pi), etc., states for that are contained 
entirely in the small space (with an analogous convention 
defining \q'), etc.). The d orthogonal small-space 
states |fc) corresponding to d selected full-space eigen- 
states |fc) are defined by 



|fc) = (p + w^) |fc) 



with 



(g|c.b) = ^(g|fc)(kb) . 



and 



M = P + Juj = P (1 + u'^uj) P . 



(5) 



(6) 



(7) 



In Eq. ([6]) the (k|p) are defined as the elements of 
the inverse of the matrix with the c?^ elements {p\k). 
With these definitions the effective operator O^s in the 
dimension-c? small space that gives the same matrix ele- 
ments as (the = analogue of) the decay operator O 
in the full space, is 



Ocff = M-^P + uj^)0{P + uj)M-^ 



(8) 



The matrix elements of this operator can be written with- 
out reference to any vectors \q) or any of the eigenstates 
fc) beyond the d that are mapped, as 



{p\ - 



(9) 



^ {p\ M-h \p,) (pi|k) {k\ O \k') (k'b2) (P2I M-^ \p') 

where the elements of M can be written in the same 
fashion as 



{p\ M \p') 



d-l 

E 



(P|k) (kb') 



(10) 



We then simply use isospin algebra to obtain the matrix 
elements of the real = 2 effective decay operator, for 
which p represents two protons and p' two neutrons. To 
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orthonormal set {\k)} that minimizes the quantity [23 
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FIG. 1. (Color online) The curves C(r), the integral of which 
gives the matrix element for neutrinoless double-beta decay. 
The solid (black) lines are the results of the full (Nmax ~ 
6) calculations with the SRG-evolved N'^LO potential, the 
dashed (red) lines are the results of the p-shell calculation 
with the effective two-body Hamiltonian and the bare decay 
operator, and the dot-dashed (blue) lines are the results with 
the effective Hamiltonian and the effective decay operator. 
The top panel is for the decay ^Ife — >■ '^Be, the middle panel 
for *He *Be, and the bottom panel for ^°He ^"Be. 



get the effective interaction we carry out a similar proce- 
dure but also include the T = states in ^Li to obtain a 
complete set of matrix elements. 

The Lee-Suzuki mapping actually is only one of many 
that are possible. As noted above, the Lee-Suzuki proce- 
dure makes the small-space eigenvectors as close as pos- 
sible to projections of the full eigenvectors without sac- 
rificing orthogonality. In other words, it constructs the 



fe=0 



lfc) - 1^) 



(11) 



This prescription seems particularly appropriate for 
a comprehensive description of the spectrum, but in 
double-beta decay we are not equally interested in all 
states. One alternative, known as the Contractor Renor- 
malization (CORE) mapping [23| is to make the small- 
space ground state |0) proportional P|0), the projection 
of the full ground state, and then construct the other 
|fc) from the set P\k) through Graham-Schmidt orthog- 
onalization. But in fact any unitary transformation of 
the |fc)'s generated by the Lee-Suzuki procedure defines 
a valid mapping. In our case, we can generate an arbi- 
trary time-reversal-preserving transformation by rotating 
the two small-space 0+ states in ^ = 6 by an arbitrary 
angle a and the two 2+ states by another angle /3. We 
will try to see whether there are values of these angles 
that are particularly suited for double-beta decay. 



III. RESULTS 

We now test the performance of our Lee-Suzuki effec- 
tive operator in heavier nuclei. Figure [1] presents our 
results for the decays '^^^'lOjjg^T.s.iOgg -^^^-^gj^ -^^g ^gg 

SRG-evolved chiral N'^LO ([IBl) interaction in a dhuj full 
space. The black (solid) curves in each of the panels de- 
note the full Nmax = 6 01/(3/3 distributions C(r). These 
curves are what the effective operators are supposed to 
reproduce. The red (dashed) curves denote the result ob- 
tained with the bare 01/(3/3 operator in the p shell, with 
wave functions produced by the effective p-shell interac- 
tion, which in turn comes from the Lee-Suzuki procedure 
for A ^ 5 and 6 discussed above. The blue (dot-dashed) 
curves are the results with the effective operator, used 
in conjunction with the wave functions from the same 
effective interaction. 

The use of the effective decay operator clearly improves 
the agreement between the p-shell C(r) and the full one 
in all three panels. One problem, however, is that C(r) is 
not itself measurable; its integral is what we want. And 
it turns out that oscillations can make apparent poor 
agreement between curves much better in the integral, 
and good agreement worse. Table |T] compares the matrix 
elements themselves for the three decays represented by 
the figure. 



TABLE I. Matix elements M fi produced by the distributions 
C(r-) in Fig.[Tl 

7 8 10 

full 1.76 0.48 0.79 
bare 1.49 0.18 0.91 
effective 1.90 0.59 1.23 
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FIG. 2. (Color online) Same as Fig. [U but for for '^He *Be 
and '^"He — ^ ^"Be only and with the CD-Bonn interaction 
conditioned for Nmax = 8 in A = 8 and Nmax = 6 in yl = 10. 



The effective operator produces a clear improvement 
in the integrated matrix element in A= 7 and (partic- 
ularly) 8, but by A= 10 the bare operator does pretty 
well and the effective operator not as well. The reason is 
apparent from the bottom panel of Fig. [1] the effective- 
operator curve, while a better approximation than the 
bare curve, is not as good when integrated because be- 
cause it is above the full curve until about r ~ 4 fm. 
The bare curve strays from the full curve at both the 
peak and dip but in opposite directions; it thus supplies 
a good approximation when integrated. 

Is this behavior a fluke? Does it depend on the shell- 
model interaction or the size of the full model space? To 
address the questions we repeated our calculations with 
different interactions and model spaces. In ^ = 7, the 
effective operator is always a decided improvement but 
in A = 8 and 10 the results are more ambiguous. Figure 
[5] and Tab, nil present results of calculations in a N^ax = 
8 space with the CD-Bonn interaction, conditioned as 
described in section |TT] for A = 8 and 10. (We do not 
show A = 7, and the size of the problem in ^°Be limits 
us to N,nax = 6 in A = 10.) Once again the effective 
operator appears to be an improvement in both cases, 
but now, as Tab. shows, the effective-operator curve 
for A = 8 cancels itself too much in the integral. And in 
^ = 10 the effective operator does better than the figure 
indicates it should. 

One might expect the procedure to work better when 
the full model space is smaller, and/or when the full re- 



TABLE II. Matix elements A4 fi produced by the distributions 
C(r-) in Fig.H 

8 10 

full -0.41 -0.67 
bare -0.48 -0.80 
effective -0.03 -0.68 

suits differ less from the bare p-shell prediction but nei- 
ther of those are entirely the case either. Figure [3] and 
Tab, mil show results of calculations in which we use the 
SRG-evolved interaction in a severely truncated space: 
Nmaa- = 2. Becausc we are using the same effective in- 
teraction as before, we have reduced the amounts of the 
wave functions that lie outside the p shell. Now all curves 
are close to each other everywhere, but the effective op- 
erator, as expected, is an improvement in A = 7 and 8. 
In A = 10, however, the bare operator is nearly perfect 
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FIG. 3. (Color online) Same as Fig. \T\ but carried out in a 
small model space with Nmax = 2. 
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FIG. 4. (Color online) C(r) for *He — >-^Be: the solid (black), dashed (red) and dot-dashed (blue) lines are as in Fig. [T] The 
dotted (maroon) and the dot-dot-dashed (green) curves result from using the normalized p— shell projections of the full ground 
states with the effective and the bare decay operator respectively. The dot-dashed-dashed (orange) curve is the result with 
bare decay operator and the Cohen-Kurath (CK) interaction [25| |. 



and the effective operator worse (though now, as Tab. If III 
shows, it is accidentally better in the integral). 

It seems, then, that except in this last example the 
effective decay operator is a decided improvement, but 
also that the oscillations in C(r) can sometimes negate 
its advantages. The natural way to do better is by adding 
three-body terms to the effective interaction and decay 
operator. We will not do so in this preliminary study, but 
can look into whether the discrepancy at the two-body 
level is due mainly to the defects in the decay operator 
or in the interaction. If we were to carry out the Lee- 
Suzuki procedure in ^ = 8 and 10, the p-shell ground 
states would resemble the normalized p-shell projections 
of the full ground states. We can therefore use these nor- 
malized projections as proxies for the states that would 
be produced by the complete A-body Lee-Suzuki p-shell 
Hamiltonian. Figure [4] shows the resulting C{t) with 
both the bare and two-body effective operators, along- 
side the curves already displayed in Fig. [TJ The perfor- 

TABLE III. Matix elements Mfi produced by the distribu- 
tions C(r) in Fig. El 

7 8 10 

full -1.06 0.70 1.10 
bare -0.90 0.37 0.99 
effective -0.92 0.45 1.14 



mance of the bare operator improves noticeably, so that 
it is about as good as the two-body effective operator 
in conjunction with the two-body effective Hamiltonian; 
its integral is = 0.63. And interestingly, perfor- 

mance gets worse when the effective two-body operator 
is used instead {Mfi = 0.74). The first result indicates 
that three-and-more-body terms in the Hamiltonian af- 
fect the matrix element, and the second that such terms 
in the decay operator do as well. 

Figure [4] has one other curve, produced by the bare op- 
erator in conjunction with the phenomenological Cohen- 
Kurath potential [2^. We can't really expect this em- 
pirical potential to reproduce the results of an NCSM 
calculation in a nucleus that isn't in reality even stable, 
but it is nevertheless interesting to see how different its 
results are. 

As mentioned in the introduction, the Lee-Suzuki map- 
ping is only one of an infinite number of possible map- 
pings. We tested for better-performing mappings in 
A = 6 by rotating the Lee-Suzuki p-shell basis vectors 
by arbitrary real angles. Though we didn't exhaustively 
explore all such possibilities, we found it difficult to im- 
prove on the Lee-Suzuki results. With the SRG interac- 
tion in Nmax = 6, rotation of the two 2+ vectors by 15° 
improved the effective operator C(r)'s marginally, but 
nowhere did we see a dramatic improvement. 
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IV. DISCUSSION 

We have shown that a nonperturbative renormahza- 
tion of the effective decay operator improves the shell- 
model's ability to reproduce Ov^P matrix elements from 
ab initio calculations. What do these results imply for the 
heavier nuclei we really care about? There is some reason 
to hope that a two-body effective operator will perform 
even better than it did here. Both QRPA and shell-model 
calculations in these nuclei show that C{r) with the bare 
operator nearly vanishes beyond about 3 fm; there are no 
oscillations of consequence in those curves [2g| . For that 
reason, the performance of the effective operator may not 
be degraded by the cancellations that play such a large 
role here. In any event it is clearly worth determining 
the effective two-body operator. 

There are several possible routes to that end. Ref. 
[27| reports coupled-cluster calculations in ^He, the first 
in a nucleus with two valence nucleons. Coupled-cluster 
techniques scale well to intermediate-mass nuclei and it 
should be possible to treat the nuclei ^^Ni, ^®Cu, and 
^^Zn with two nucleons in the /5/2TO9/2 shell (along with 
the simpler nuclei with one nucleon in that shell). Such 
calculations will play the same role as the NCSM calcula- 
tions here, allowing the extraction of a two-body effective 
Hamiltonian and decay operator, which can then be used 



in a shell-model calculation of the decay of e.g. Ge. If 
the Hamiltonian proves inadequate, a better phenomeno- 
logical one can be substituted, though at the obvious 
cost of consistency. Other techniques, including many- 
body perturbation theory with soft renormalization- 
group- produced interactions l28i and the in- medium sim- 
ilarity renormalization group |29| offer a path to effective 
operators as well. 

For a truly accurate calculation we probably will need 
three-body decay operators. All the methods mentioned 
above can produce these in principle, but at least a few 
years of development are required. In the meantime, we 
can examine the question of whether two- and three-body 
decay operators will be sufficient (or whether horribly 
complicated four-body operators will be required) by ex- 
tending the p-shell tests we report here. Work in that 
direction is underway. 
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